> N85-31199 


MATRIX TRANSFER FUNCTION DESIGN FOR 
FLEXIBLE STRUCTURES— AN APPLICATION 

T. J. Brr— mi, A. V. Ciwnpito, A. L. Dam, C. L. Cimfrfiw, and C. L. Wong 

A -nspace Corporation 
Lite. /%nge>es, CA 90009 


ABSTRACT 

The application of matrix transfer function design techniques to the 
problem of disturbance rejection on a flexible space structure is demonstrated. 
The design approach is based on parameterizing a class of stabilising 
compensators for the plant and formulating the design specifications as a 
constrained minimisation problem in terms of these parameters. The solution 
yields a matrix transfer function representation of the compensator. A state 
space realization of the compensator is constructed to investigate performance 
and stability on the nominal and perturbed models. The application is made to 
the ACOSS (Active Control of Space Structures) optical structure. 


I. INTRODUCTION 


The problem of flexible space structure control has motivated a great 
deal of research for theoreticians and practitioners of multivariable control 
design. In spite of the efforts directed in this area there still remains a 
significant gap between the multivariable theory and the control design 
implementation. This gap stems from two sources. The first difficulty is one 
of problem specification. Translation of complex system* requirements and 
constraints into the specific mathesmtical cost functionals required by most 
design methods may be impossible in many cases. Free parameters in the chosen 
design methodology may not be traceable to the parameters which describe the 
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system in terms of desired performance, plant uncertainty, hardware 
limitations, etc* A second roadblock to the implementation of modern control 
design techniques is the iack of reliable algorithms and software to perform 
the sophisticated mathematical manipulations required by these techniques. 
Recent years have shown very considerable advances in this field (see [1]) but 
much remains to be done. 

Most of the MIMO (multi- input/multi-output) compensators which have 
actually left the textbook and been calculated in computers are based on state 
space methods, and, in particular, LQG ( Linear-Quadrat ic-Guas a ian) design 
theory. This is due in part to the long history of development of these 
design techniques as well as the availability of reliable algorithms to solve 
matrix Riccati equations and the ease of performing most state space 
manipulations. Frequency domain techniques for calculating MIMO feedback 
systems have been avoided. The extensions of classical frequency domain 
concepts to MIMO systems have not been totally satisfying and calculations 
involving matrices of transfer functions present an entirely new set of 
problems. Nonetheless, frequency domain design is still appealing and certain 
feedback notions cannot be adequately expressed without reference to transfer 
functions. 

We have carried through a compensator design for a flexible structure 
based on transfer function parameterization techniques. General theories of 
feedback control system parameterization have been developed by several 
authors ([2], [3], and [4]). The goal of a parametric approach is tne 
selection of a set of numerical quantities, along with an acceptable range of 
values, which span a class of possibly acceptable compensators and, wi - 
which, one is able to adequately express the system requirements in terms of 
cos 1 8 and constraints. A particularly simple parameterization for stable 
plants was introduced by Zanies, [4], and exploited for the unity feedback 
conf iguration of Figure l bv Desoer and Chen [5]. This is the 
parameterization we will implement here. The details are in section IV. 
Previous examples of this design approach can be found in [6] and [7]. 
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II. ACOSS STRUCTURE 


The ACOSS optical structure was developed by the Charles Stark Draper 
Laboratories as a control design test specimen to evaluate the design 
approaches developed for the DARPA ACOSS program, f 8] • It was designed to 
exhibit the closely spaced, low frequency mode distribution expected on some 
future space systems. The structure is provided as a finite element model 
having 84 dynamic degrees of freedom (see Figure 2). In addition to the 
nominal structure, two perturbed structures were defined to represent plant 
uncertainty. The perturbed models represent mass and stiffness variations of 
approximately 10Z. The nominal model is denoted PO, the perturbed models are 
P2 and P4. 

The performance goal is expressed in terms of a line of sight error on a 
focal plane on the lower section of the truss as shown in Figure 2* The error 
has two angular components and a defocus component resulting from deviations 
in the optical path due to structural vibrations. Three rigid mirrors 
determine the optical path. These are assumed to be rigidly mounted to the 
structure. Two disturbances are defined ou he structure as shown in the 
figure. For our design problem we are only considering the disturbance 
propagating from the equipment panel and we assume it has a flat PSD out to 
5 Hz. The equipment panel is isolated from the structure by a spring-damper 
system. The residual disturbance propagation through this isolation system 
into the line of sight is still unacceptably high. The control problem is to 
further reduce this residual with active structural control. 


III. MODEL SELECTION AND ACTUATOR PLACEMENT 


For the current design problem we chose a 5 mode model of the structure, 
selecting those modes having most significant influence between disturbance 
and line of sight. The modal influence was determined based on ideas from 
internally balanced coordinates. For a description of internally balanced 
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coordinates see f 9 ] and for an application to modal coordinates see [10l. 
Given a second order model description. 
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Using the modal disturbance influence matrix for the g^'a and the line 
◦f sight measurement matrix for the h^'s the 5 highest rank modes are 
tabulated in Table 1. Agreeing with our intuition, these turn out to be two 
isolator rotations, two isolator translations, and the first bending mode of 
the upper truss* A description of the modes of the structure can be found in 
[ 8 ]. 
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Table 1* Design Modes 


The line of sight measurement matrix is a function of 21 nodal degrees 
of freedom, from among these 21 degrees of freedom we chose to locate three 
force actuators (assumed to be of the momentum exchange or proof mass type) to 
control the three line of sight measurements* To make this selection an 
appeal is again made to the approximate second order modes of equation (3). 
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If the forcing function on the right of (l) is g .u where g.. is the 
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influence of the j actuator, j * 1,...,21 on the i mode, i * l,...,5, 
then we denote the corresponding second order mode by and define 
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Here, a. is a measure of the influence of the j** 1 actuator on the 
line of sight for the selected 5 mode model. We chose three actuators whose 
force directions span the three spatial directions and have large a. with 
respect to the total 21 possible actuators.* Two of the actuators selected are 
located on the corners of the primary mirror and the third is on the lower 
truss • 


To complete the description of the design plant we assumed the 
availability of direct measurements of the line of sight. No other sensors 
were used for the control design. We now have a state space description of 
the design plant in modal coordinates, 

x - Fx + Gu + Dd (5) 

y * Hx (6) 


where u is .he actuator command and d is the disturbance input. 

Per calculation of the compensator we need a transfer function 
represents ion of the design plant. The convenient representation for 
constructs g state space realisations of the compensator is a polynomial 
matrix coprime factorisation [11, 12 1 , that is, P * ND ^ where N and D are 
coprime polynomial matrices. An algorithm to construct a coprime 
factorisation from a state space description can be found in [13]. 
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IV. DESIGN PROBLEM 


The feedback configuration used for the design is shown in Figure 1. 

The closed loop system is referred to as P is the open loop design 

plant , a 3x3 transfer function given from the state space equations by 

H(sl-F) *G. The inputs are u^ and u^ with the reference input, u^, 

identically zero. The outputs are and y ^ with the line of sight 

represented by y . The disturbance propagates into the line of sight 
1 ~~ , -I 

through the transfer function P =* H(sl-F) D and may thus be represented as 

*v 

an additive disturbance, Pd, at the plant output. 

The closed loop system transfer function is defined to be 


H 


yu 



(7) 


Stability of H can be taken to be closed loop stability. 
yu 

expressed in a simple parameterized form as 


H may be 
yu 


Q ~QP 

H 

yu 

PQ P(I-QP) 


( 8 ) 


where Q is referred to as the Zames parameterization, [4], with 

Q * C(I + PC)" 1 . 


(9) 


We state here the fundamental result from [5] which is the basis of this design 
approach. 

Fact : For P exponentially stable and strictly proper, Q is exponentially 
stable and proper if and only if 
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(i) C is proper and 

(ii) Hy u is exponentially stable and proper. 

When this is the case the compensator is given by 

C - Q(I - PQ)" 1 (10) 

In other words designing stabilizing compensators for is equivalent to 
specifying exponentially stable, proper Q. 

From (8) we see that the I/O map, that- is, the transfer function from 
U 1 t0 y 2 * s 


H * H = PQ (ll) 

y 2 U l 

Given an invertible plant transfer function, P, one can see from the 
reLation (11) that a parametrization of the closed loop system by Q is 
equivalent to a parameterization by H. Moreover, for P exponentially stable, 

Q exponentially stable implies the same for H. But since 

Q - P~ l H (12) 

it becomes clear that exponential stability of H only implies exponential 
stability of Q when P has no unstable zeros. However, by imposing an 
additional condition on H, namely that H has the same right half plane zero 
structure as P, then parameterization by such H is equivalent to parameter iza- 
tion by exponentially stable Q. If a proper compensator is desired the 
additional constraint of properness of Q is required and will result in an 
excess pole over zero constraint on H which depends on P. 

Parameterization by the I/O map, H, may simplify the design problem and 
allow the designer to more directly specify his design objective. For 
example, for a disturbance attenuation problem, the closed loop disturbance to 
output map, or sensitivity map, is simply given as (I - H). In addition, in 
some applications, a decoupled I/O map is desirable and one is directly able 
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to parameterize a diagonal H. This is the approach we take for this design. 
Calculating the transmission zeros of cur design plant using the QZ algorithm 
[14] we find that there are no zeros in the right half plane so we may freely 
specify H as diag(h^, i=l,2,3) with each h^ of the form 
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where g is a gain and 
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This parameterization has 21 parameters consisting of the gains , and 
second order damping and frequency terms* 

We set for ourselves a design goal of minimizing closed loop response to 
the disturbance over a low frequency band of 5 Hz* To achieve this we define 
a constrained optimization problem as follows: 

Minimize 


subject to 
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:Low frequency rcise rejection 
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The matrix under the norm of J is diagonal so we aimj^v take the 
Euclidean vector norm of the diagonal. The minimization o* ne cost at 5 Hz 
and the DC unity gain constraint will result in disturbance rejection across 
the 5 Hz band* The P term in the cost weights the diagonal terms in (I - H) 
according to the way the disturbance propagates through the structure. 

In general (I - H) is the ratio of the relative uncertainty in the I/O 
map to the relative uncertainty in the open loop plant. More nreciseiv 

(AH)H* 1 - (I~H)(AP)P” 1 (17) 


where 


AP » P - P 

(ltJ7 

AH - H - H 

(19) 


for a "perturbed" plant P which results in a perturbed I/O map H. In effect, 
minimization of J reduces the impact of plant uncertainties on closed loop 
system performance. 


Having specified the optimization problem one can use numerical or 
analytical means to solve it. Omitting the details, we calculated a local 
minimum to this problem analytically. The achievable performance is clearly 
dependent on the bandwidth, 0)^, For a 
satisfies 
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Given this solution we can adjust the bandwidth of each of the three 
loops to achieve a desired performance level. To achieve 0.04 reduction in 
each channel we have the following parameter values: 
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V. COMPENSAiOR REALIZATION 


Having arrived at parameter values we have specified tne desired I/O 
map, H. The compensator which will produce this l/° mnp is 


C - Q(I - PQ) 


(24) 


* p“Hi(l - H)" 1 . 


(25) 


Since the I/O map is given by diag(h^ ) with each h^ of the form 
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( 26 ) 
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the compensator becomes 


C 


P *diag 



(27) 


We have already expressed P as a polynomial matrix coprime factorization, 
P - ND Thus (27) becomes 


C * 


D N l diag 



(28) 


Since the degree of d^ - n^ is 4, we can factor this polynomial into two 
quadratics as 


d. - n. « d. d. , i - 1,2,3. 

ii i L i 2 


Hence (29) can be rewritten as 


C 


D { diag(d. )N } * 



(29) 


(30) 


By inspection, diag(d. )N is column- reduced [12], and has column degrees 

H - . 

equaling those of D. Consequently D {diag(d. )H} is proper and has a state 

l l 

space realization [12, Sec. 6.4]. Now, since diag(n./d. ) also has a state 

1 V 

space realization, the two realizations can be cascaded to yield a realization 
for C. 


57 



VI. RESULTS 


Having computed a state space description of tho compensator we are now 
able to determine closed loop stability for various versions of the plant 
simply by extracting the eigenvalues from the closed loop state equations 
derived from Figure 1. We find that for all three versions of the plant (PO, 
P2, and P4) the five mode description remains stable under feedback by our 
compensator. 


To investigate the robustness of the design with respect to unmodeled 
dynamics we appended additional modes to the plant model and found that the 
closed loop system became unstable in almost all cases. Upon investigation of 
this problem we discovered that though the 5 mode design plant was minimum 
phase, the addition of almost any other mode or set of modes resulted in a 
nonminimum phase plant. Information about these unstable zeros was not 
available in the design plant so the resulting compensator tended to place 
closed loop poles at these zeros. Thus the stability problem experienced is 
one of modeling or model reduction. In general, any control design approach 
must have information about the righ*^ h*lf plane zeros of the plant. 


The performance of the closed loop system remained very consistent with 
the predictions made during the design stage. The steady state RSS response 
at 5 hz of the two angular components of the line of sight is given as a 
fraction of open loop response for the three models by: 
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PO 

P2 

P4 

-2 

-2 

-2 

4.3 X 10 

5.1 X 10 4 

4.6 X 10 


The broadband disturbance attenuation is illustrated on the Bode plots 
of Figures 3 and 4 which compare open and closed loop response. Across a 
significant portion of th* 5 Hz band the performance improvement is 3 to 4 
orders of magnitude. 
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VII. CONCLUSIONS 


We have demonstrated the applicability of a transfer function parameteri- 
zation design approach for problems of broadband disturbance attenuation on 
flexible space structures. This methodology provides the control designer 
with a great deal of flexibility to meet system requirements by the choice of 
parameter set and selection of cost function and constraints. Although the 
implementation of this technique requires difficult numerical calculations 
involving matrix transfer functions, algorithms and software for these types 
of problems are already emerging. The success of this approach is dependent 
on an appropriate parameter selection in which to express the problem 
specifications. This suggests research, probably application specific, which 
addresses the issues of problem description and requirements interpretation in 
the control design process. 
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